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ABSTRACT 



A quasi-periodic signal is sampled at a specified rate, 
and then a predicted value of the signal is computed 
from a set of 2M+1 time lagged signal samples. The 
time lagged samples are centered in time at an integer 
multiple P of the signal's sampling period T Si where 
P-T* is approximately one period of the input signal. 
The predicted signal value is computed by multiplying 
each of the 2M-j-l time lagged samples by a corre- 
sponding predictor coefficient c(i) and then summing 
the resulting products. The predicted signal value is 
subtracted from the actual signal value to obtain an 
error signal €. During each successive sampling period, 
the predictor coefficients are updated by adjusting the 
previously computed predictor coefficients by an 
amount proportional to the error e multiplied by each of 
the 2M+ 1 time lagged signal values. Using the updated 
coefficient values, a phase delay is computed, and then 
the signal's period is computed as the sum of phase 
delay and P-Tj> which is the average integer time lag of 
the sampled signal values. When the phase delay is 
greater than one half of the sampling period T 5) the 
value of P is increased or decreased, as necessary, to 
keep the time lagged signal samples approximately one 
period delayed from the current signal sample. 

17 Claims, 3 Drawing Sheets 
^122 



Quasi-Periodic 
Input 

Signal — ■» 



120 



124 



DSP 



Sample 
&Hold 

10KHz Clock 



126 



Sample Storage Memory 



144 



132 



146 



/M40 



x(k) 



AMDF 
Routine 



Coefficient 

Update 

Routines 



SIN, COS, 

ARCTAN 

Tables 



I Multiplier I 



Iff! 



EH 

£(k)| 



Pitch, 

MIDI Commands 



130 



CPU 



/M42 



Coefficients, 
Pitch, Delay, e, 
Average Energy 



Commands , f 



Pitch, Delay, e, 
Average Energy 

if 150 



Microprocessor 



03/12/2004, EAST Version: 1.4.1 



U.S. Patent 



Oct 4, 1994 



Sheet 1 of 3 



5,353,372 



^10 2 ^104-1 ^104-2 



• • • 




104-2M 



-1 _ 



C(M) 



x(k) 



*(IK-*e(k) 

-^ -J v f11 2 



{C(i)}« 
i=-M to M 



FIGURE 1 



Coefficient 

Update 

Computation 



Quasi-Periodic 
Input 

Signal — — 



120 



122 



124 



DSP 



f^Pl 6 — I Sample Storage Memory 

& Hold ■ i ii i 
10KHz Clock 



^126 



144 
132 

146 



140 



x(k) 



AMDF 
Routine 



Coefficient 

Update 

Routines 



SIN, COS, 

ARCTAN 

Tables 



128 

Vector U^K 



Multiplier 

hhe 



•108 



+ I 



e(k) 



Pitch, 

MIDI Commands 



v ^130 



CPU 



JL 



142 



Coefficients, 
Pitch, Delay, e, 
Average Energy 



Pitch, Delay, £, 
Commands .Average Energy 



Microprocessor 



FIGURE 2 



03/12/2004, EAST version: 1.4.1 



U.S. Patent Oct. 4, 1994 sheet 2 of 3 5,353,372 



f200 

INITIALIZATION: 

Compute Initial Period Using AMDF Technique: 
Period s TO 3 P-Sampling Period) 
Set c(i)=0 for all i, except c(0)=1 . 



Compute: 



/201 



e 2 >px 2 

Indicates 
Loss of 
Pitch 
Tracking 



Error Signal, e^, Mean Square Error, e 2 



202 



Average Signal Power, x 2 



204 



C k+1 = C k + 



(2M+1)x* 



~ x k E k 



2flfi 



<p = 2n/(Previously Computed Period) 



208 



6=arctan 



f M 

Z(c(i)-c(-i)) sin(coi) 
i=1 



M 

,c(0)+ Z(c(i)-c(.i))cos((«) 
1=1 



210 



Phase Delay = 



CO 



212 



Period " Sampling Rate + P"«eDelay 
1 



Pitch 



Period 



2H 




f220 

Increment/Decrement P, 
Set c{i) = c(-i) 
Start Hold Period 



FIGURE 3 



03/12/2004, EAST version: 1.4.1 



U.S. Patent Oct 4, 1994 sheet 3 of 3 5,353,372 



Quasi-Periodic 

Input 

Signal 



252 



Pitch Tracker with P=P0 - 1 



/-254 



Pitch Tracker with P=P0 



r 256 



Pitch Tracker with P=P0 + 1 



Pitch, 

MIDI Commands 



Commands 



260 



Microprocessor 



FIGURE 4 



250 



Quasi-Periodic 

Input \ 

Signal 



272 



Pitch Tracker with P=4 



274 



Pitch Tracker with P=6 



276 



Pitch Tracker with P=8 



JL 



278 



Pitch Tracker with P=202 



Pitch, 

MIDI Commands 



Commands 



280 



Microprocessor 



270 



FIGURE 5 



03/12/2004, EAST Version: 1.4.1 



5,353,372 

1 2 

Then autocorrelation of the input signal is computed as 
ACCURATE PITCH MEASUREMENT AND follows: 
TRACKING SYSTEM AND METHOD 

q+N-\ (Eq. 2) 

The present invention relates generally to pitch de- 5 xxx{m) - ^ 2 *0Xi + m ) 
tection systems and methods for detecting the pitch of 

quasi-periodic sound sources. where q represents a starting point in the input signal, 

BACKGROUND OF THE INVENTION and N is an integer representing the number of signal 

samples to be used in the computation. The main peak in 

Pitch detection is of interest whenever a single quasi- 10 ^ autocorrelation function is at the zero-lag location, 
periodic sound source is to be studied or modeled. For m= o. The location of the next peak gives an estimate of 
instance, the trajectory of a sound's pitch, also called fa e signal's period, and the height of the autocor- 
the fundamental frequency, over a period of time can relation function given an indication of the periodicity 
also be used to synthesize similar or related sounds using Q f tne signal. Disadvantages of this method are (1) that 
various speech or music synthesis techniques. An exam- 15 it is computationally intensive, making real time opera- 
pie of a quasi-periodic sound source is a singer's voice tj on difficult or impossible, and (2) it only gives a rough 
singing a particular note (e.g., high C). The sound gen- estimate of the signal's period with an accuracy no finer 
erated by the singer typically has a certain amount of m time resolution thant the sampling period. One ad- 
vibrato, or pitch modulation, noise and aperiodichy in vantage of the autocorrelation method is that it gener- 
the wave shape, making the sound quasi-periodic rather 20 a tes a pitch value without requiring any prior informa- 
than a pure periodic signal. tion about the approximat e value of the pitch. 

Prior art pitch detection methods can be classified — . — ; ; 

into three categories: frequency domain techniques, - Average Magnitude Difference Function 

time domain techniques, and methods which use both ^ Another prior art time domain pitch tracking method 
techniques. The present invention is a time domain is called the Average Magnitude Difference Function 
detection technique. (AMDF). The AMDF pitch detection method is the 

Important characteristics of a commercially useful complement of the autocorrelation function, in that it 
pitch tracker are (1) real time operation on signals in the measures the difference between the waveform and a 
0 to 5 KHz range using inexpensive hardware, (2) accu- 3Q time lagged version of itself. The AMDF function is: 
racy, and (3) generation of pitch values at regular (uni- 
form) time intervals. Another useful characteristic is the q + N _ x (£q. 3) 
ability to detect when the pitch tracker has lost track of AMDF{m) = ^ 2 [*(o - *0* + m)| 
the input signal, as may happen when the pitch of the 

input signal suddenly jumps, and to recapture the input 35 lagj m=0) potion of the AMDF function is 

signal without human intervention. identically zero, and the next significan^nuins„a.Ukely 

In time domain "feature detection" methods, the esti mate of the period. Other nulls will occur at integer 
input signal is usually preprocessed to accentuate some multiplis^ 

time domain feature, and then the time between occur- using techniques known to those skilled in the art, to aid 
rences of that feature is calculated as the period of the 40 m detection of the first null. A primary advantage of 
signal. The pitch and period of the input signal are this pitch detection method is that has a smaller compu- 
related by the equation: tational burden than the autocorrelation method, be- 

cause it uses subtraction where autocorrelation uses 
l (Eq. l) multiplication. 

Period 45 The difficulties of the AMDF pitch detection method 

arise from the finite sampling rate of the signal, noise, 
A typical time domain feature detector is implemented and signal stationarity. If the signal is truly periodic 
by low pass filtering the signal, then detecting peaks or with period To, and To is an integer multiple of the 
zero crossings of the filtered signal. Since the time be- sampling period T St then all nulls at integer multiples of 
tween occurrences of a particular feature is used as the 50 To will be identically zero. However, if the period of the 
period estimate, feature detection schemes usually do signal is not an integer multiple of T s , the first null 
not use all of the data available. Selection of a different (m^=0) actually exists between two values of m. For 
feature often yields a different set of pitch estimates. example, if the sampling period is 0. 1 milliseconds, and 
Since estimates of the period are often defined at the the input signal has a pitch of 440 Hz and a period of 
instant when the features are detected, the frequency 55 0.002272727 seconds, then the true first null will be 
samples yielded are not uniformly distributed in time. between samples m=22 and m=23. A coarse estimate 
To avoid the problem of non-uniform time sampling, a of pitch is tolerable for many speech applications, but is 
window of fixed size is moved through the signal, and a not acceptable for analysis and synthesis of music. Corn- 
number of detected periods within each window are pared to the small computational burden of computing 
averaged to obtain the period estimate. 60 the AMDF, there is no known economical method of 

accurately interpolating between samples to find the 
true period. This implies that the sampling rate must be 
Other prior art time domain methods include the use sufficiently high to yield the high accuracy required for 
of autocorrelation functions or difference norms to musical applications. 

detect the similarity between the waveform and a time 65 If the signal is quasi-periodic (e.g., amplitude modu- 
lagged version of itself. The autocorrelation method lated, frequency modulated, or corrupted by noise), the 
works as follows. The input signal is sampled and stored nulls will never be zero, even if To is an integer multiple 
at discrete time intervals, such as every 0. 1 milliseconds. of T s . The problem of interpolation between lag samples 
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Autocorrelation Function 
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to obtain an accurate piteh estimate is even further DESCRIPTION OF THE PREFERRED 

complicated in the case of a frequency modulated sig- EMBODIMENT 
nal. 

The present invention is an improvement on the Referring to FIG. 1, the pitch tracker of the present 
AMDF pitch detection method. Using a FIR (finite 5 invention is conceptually a tapped delay line imple- 
impulse response) filter model, the present invention mented using a sequence of delay elements 102, 104 for 
provides a methodology for computing a signal's period generating time lagged, representations of an input sig- 
with great accuracy. nal x(k). The first delay element 102 has a delay length 

of P— M where P is an integer number of time sample 
SUMMARY OF THE INVENTION 10 periods representing a rough estimate of the signal's 

In summary, the present invention is a system and P eriod - For example, if the input signal were a 440 Hz 
method for tracking the pitch of quasi-periodic signals. tone » and P eriod T * were ai Seconds P 

The quasi-periodic signal is sampled at a specified rate, w °? d be ™> presenting an approximate period of 2.3 
and then a predicted value of the signal is computed 15 ^iseconds. The number of subsequent umtory time 
from a set of 2M+1 time lagged signal samples. The dda / s ™ » where M is usually a .^integer 
. 4 . . * j • i: * * * such as 1, 2 or 3 and preferably less than eight. It should 

Ume lagged samples are centered m time at an integer be ^ £ e mvention can ^ med ^ ^ 

multiple P of the signal s samphng period T, where f vdue ^ ^ values of M 

PT, is approximately one penod of tfie input signal. ^ J^tfonal computational expense. 
The predicted signal value is computed by multiplying 20 A sum of 2M + 1 samples of the time lagged 

each of the 2M+ 1 time lagged samples by a cane- signa]j at ^ ddays ranging from P-M to p+ M> ^ 
spending predictor coefficient c(i), and then the pre- formed multipliers 106 and summer 108. The 

dieted signal value is subtracted from the actual signal we i g hts applied to each time lagged samples are called 
value-to obtain an error signal e. predictor coefficients c(-M) to c(+M). The resulting 

During each successive sampling period, the predic- estimated signal value is herein called x(k): 
tor coefficients are updated by adjusting the previously 

computed predictor coefficients by an amount propor- M (Eq. 4) 

tional to the error € multiplied by each of the 2M+ 1 *(*) = M x(k - P + iMO 

time lagged signal values: 30 i== 

c =c +2uX*m As will be discussed below, the predictor coefficients 

* + i k + yjLM are selected so that the estimated signal x(k) matches the 

where C*is a vector containing the 2M+ 1 predictor 38 ^ly as possible. 

" K . - j i - . j The estimated signal value x(k) is subtracted from the 

coefficients computed for toe penod k, 2^ is an adap- 35 

tion parameter that controls the rate at which the pre- ^ ^ m HG ^ fe _ 

diction coefficients are adjusted and X* is a vector ^^^ ooa ^^ m 

contaimng the 2M+ 1 time lagged signal values at time ^ m to date ±e yalues of ^ predictor co^. 
period k. Using the updated coefficient values, a phase ents - for j= _ M tQ +M Methods for computing 
delay is computed, and then the signal's penod is com- 40 and updating these coefficients will be discussed below 
puted as the sum of phase delay and P-T,: m ^ section entitled "Computation of Predictor Coef- 

ficients". For now, it should be assumed that the pitch 
tracking system computes the optimal coefficient values 

, , , , M& , ir A ^ t t a so as to ininimize the mean square error (MSE)?(also 

When the phase delay is greater than one half of the 45 „ , - . , ^ x , . , «„JL™ 

v <i +v. & -r\ • * _ , _ . called the error signal power), which is the average 

samphng penod T St the value of P is increased or de- . * r *T ^ 7 . , - x 

k & *~ * , ,. , * i value of the square of the error signal €(k). 

creased, as necessary, to keep the tune lagged signal ^ 

samples approximately one period delayed from the Digital Signal Processor Implementation 

current signal sample. ^ Referring to FIG. 2, the preferred embodiment 120 of 

BRIEF DESCRIPTION OF THE DRAWINGS the pitch tracker of FIG. 1 uses a digital signal proces- 

. . . .„ sor 122, such as the 56001 made by Motorola. The DSP 

Additional objects and features of the mvention will 122 fr^es ^ m ^ og to digital converter (ADC) 124 

be more readily apparent from the following detailed for sampling a qua si-periodic input signal at discrete 

description and appended claims when taken in con- ^ ^ mtervalS( such ^ once every 0.1 inilliseconds (i.e., 

junction with the drawings, in which: at a rate of j 0 siglia i samples are stored by 

FIG. 1 is a block diagram of an FIR signal predictor ^e DSP in signal storage memory 126, which can either 

in accordance with the present invention. ^ e a sm ft register or a block of random access memory 

FIG. 2 is a block diagram of the hardware of a pre- configured as a circular buffer, or a combination of 

ferred embodiment of a pitch tracking system in accor- ^ ^e^ two (i.e., a shift register can be used for the last 

dance with the present invention. 2M+1 elements of the memory 126, while random 

FIG. 3 is a flow chart of a preferred embodiment of access memory is used for the first P— M elements of 

the pitch tracking method of the present invention. the memory 126). The size of the signal memory 126 is 

FIG. 4 is a block diagram of a pitch tracker using governed by the lowest frequency signal that may need 

three parallel pitch tracking subsystems. 65 to be detected. For instance, if input signals with pitch 

FIG. 5 is a block diagram of a pitch tracker system as low as 20 Hz may be received, and a 10 KHz sample 

using a large number of parallel pitch tracking subsys- clock is being used, the signal storage memory would 

terns. need to be able to hold approximately 510 signal sam- 



Period=P-T,+ Phase Delay 
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pies (i.e., 500 samples for one period, plus M additional 

samples for use by the FIR filter). ^ (Eq. 7) 

Next, a vector multiplier 128.is used to multiply the y ^ 

last 2M+1 time lagged signal samples by their corre- 
sponding predictor coefficients, and a summer is used to 5 By adding the computed Phase Delay to the time delay 
subtract the resulting 2M+ 1 values from the current P* T * associated with the center delay element of the 
signal sample value x(k) to produce an error value c(k). filter, the net time delay of the predictor is computed. 
During each sampling period, a central processing unit This total delay is then used to compute a period and 
130 in the DSP receives the error value c(k) as well as 1Q frequency estimate: 
a vector X* comprising the 2M+ 1 time lagged signal 

values, and then computes a new set of predictor coeffi- Period =T 0 -P.T jr +Pha Se May (Eq, 8) 

cients using one of a number of coefficient update rou- Frequency=F 0 = i/T 0 (Eq. 9) 

tines 132 stored in the DSP's software memory 140. The 

resulting coefficient values, as well as a computed per- 15 where Xjis ^ period in seconds and T 0 is the 

iod, pitch, and average signal energy value are stored in period estimate. 

a parameter memory 142. Since Tq ^ not prec isely known when the pitch 

In the preferred embodiment, when the DSP 122 is tracker first begins its initial calculation, the initial value 
initialized upon power up or upon being reset, it com- 0 f o is somewhat coarse, which in turn causes the first 
putes an initial approximate Period value equal to an 20 comput ed value of the FIR filter phase 0 to be some- 
integer P number of sample periods using the AMDF what ^,3^ jf a more accurate estimate is required, the 
technique described above. Thus, an AMDF routine last com puted Period value is used to compute o>, or the 
144 stored in the software memory 140 is executed by calculation of a> and the pitch estimate is iterated until a 
the CPU 130 at initialization, thereby generating an 25 de sired accuracy is reached. 

initial period value P. The number of computations required to compute the 

Also stored in the software memory 140 are sine, FIR filter phase 9 can be reduced by exploiting the 

cosine and arctan tables 146 which are used in certain symmetry of the cosine and sine functions, as well as the 

computations, described below, to calculate sine, cosine symmetry of the filter definition. Equation 5 reduces to: 

and arctan values by interpolated table lookup. 30 

Finally, for reasons that are described below, a mi- 
croprocessor 150, such as the 68000 made by Motorola, ^ 

is used to monitor the DSP's operation by monitoring 

parameters such as the computed pitch, phase delay, 

integer delay P, average signal energy, and error signal 35 © = arctan 

e(k). As will be described below, when the monitoring 

microprocessor 150 determines that the absolute value 

of the phase delay in the FIR filter exceeds half of a _ ^ , , 

sampling period, it commands the DSP's CPU 130 to For further computational savings, sine, cosine and 
either increment or decrement the integer delay P so as 40 ™ vaIues can be calculated by mterpolated table 
to keep the FIR filter centered on the signal's period. lookup. 

The monitoring microprocessor 150 also can be used to Computation of Predictor Coefficients 

periodically generate MIDI commands (to be sent to a . . . .. « , . , ^ „. ,„ . 

.« j * 1 ct +1 ♦ii The key to keepmg the pitch detector tuned to the 

musical synthesizer) or other data signals reflecting the A . . . . jf ^ , ^ , 

pitch value detected by the DSP 122. 45 ™ d to accurately computmg the penod 

p y of that signal is computmg a set of predictor coefficients 

Determining the Phase Delay of the FIR Filter that minimize the mean square error (MSE)€ 2 . The 

The phase 9 of the FIR filter, relative to the Pth * a Potation of three different technique 

delayed sample, is a function of the coefficients: «, for *ese coefficents, herein called I the Co- 

J r variance method, the Recursive Least Squares (RLS) 

method, and the Least Mean Squares (LMS) method. In 
the preferred embodiment, the LMS method is used 
because it is reasonably accurate and because it requires 
55 the least number of computations, which enables real 
time performance using modest equipment. 

The Covariance method of computing prediction 
coefficients works as follows. The coefficients c(i) are 
The angular frequency o> in the above equation is the assumed to implement a least squares prediction of the 
frequency of the input signal. In other words, ^ sign's value at a point P steps ahead of a selected 

point in time. A set of coefficients is computed over a 
= 2n _ 2v (Eq. 6) defined data set, such as a set of 128 signal samples. The 

To ~ PTs coefficients solve the following linear set of equations: 

where To is the actual period of the input signal and $5 RC=r (Eq. 11) 

P-Tj is a rough estimate of the input signal's period. 

Using the above equations, the phase delay of the where C is the vector of predictor coefficients, r is the 
filter is computed as: P-delayed covariance vector, 



? MO " e(-/))sin(a») 
1=1 

c(0) + f (c(0 - <t-i))cas(oi) 



0 = arctan 
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. X c(i)sin(o)/) 
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+ p + 1 - mi 



£[x(i>(r + P + M)) 



and R is a square matrix of expectation values: 



8 



inversion can be updated recursively using the follow- 
ing matrix inversion lemma: 



to get the following equation: 



R k = — r 



(Eq. 15) 



(Eq. 16) 



10 - MX + aJr^li^fc) 

where X=l— a. With this estimate of the autocovari- 



R = 



£[*0X0J +1)1 - • -mOxV + 2^)3 

£W + Wffi £[*(' + Wf +1)] • • •£[*(' + + 2Af)] 



£!*(/ + 2Af)x(i)] + 2A0x(i + !)]•• E[x(i + 1M\x{i + 2J0] 



where the expectation function EQ here is defined as: 



ance matrix, the predictor coefficients are updated 
25 adaptively: 



EM + a W + A)] = jj *X J *</ + o)xti + *) 



(Eq- 12) 



If the matrix R is invertible, the coefficient vector C is 30 
given by 



c=R- 



(Eq. 13) 



35 



The input signal is processed in blocks of length N 
(e.g., 128 samples), yielding a set of predictor coeffici- 
ents and thus a period estimate, for each block. Calcu- 
lating each set of coefficients requires 2NM multiplies 
to compute the covariance matrix and on the order of 
M 3 operations for the matrix inversion. Since M is small, 
the size of the R matrix is small and therefore most of 
the computational burden is associated with computa- 
tion of the elements of the covariance matrix R. 

The Recursive Least Squares coefficient computation 
method implements a recursive approximate Newton's 
method solution of equation 13. The covariance matrix 
R is approximated by an exponentially decaying sample 



R k =(l-a)R k -i+aX k X k T 



(Eq. 14) 



50 



where a is a tunable decay parameter with a value 
between 0 and 1, 



**+ 1) 



and 



65 



The covariance matrix at this point can be inverted 
directly. To reduce computation, however, the matrix 



C k +x~C k +aR-hXk<k) 



(Eq. 17) 



The number of operations required to perform the 
RLS method is on the order of M 2 , but since M is typi- 
cally low, other computational factors^ usually domi- 
nate. To avoid numerical problems, the R- 1 matrix may 
have to be updated using matrix factorization tech- 
niques. 

The Least Mean Squares method is a gradient steep- 
est descent methodology using the instantaneous error 
€(k) to estimate the gradient of the error surface. Each 
coefficient is adjusted each sample period by an amount 
proportional to the product of (A) the instantaneous 
error e(k) and (B) the signal value which is associated 
with the coefficient being adjusted. This corresponds to 
setting the R~ 1 matrix to the identity matrix in the RLS 
coefficient equation 17, and yields the LMS update 
equation: 



45 



(Eq. 18) 



The adaption parameter 2u. replaces a by convention 
from the Newton's method derivative (the factor of two 
comes from taking the multi-dimensional derivative of 
the error function), and controls the dynamics and sta- 
bility of the process of adapting the coefficient values. 
Stability is ensured if: 



55 



0<ji<((23/+l)x2)-i 



(Eq. 19) 



where x 2 is the signal power. The adaption parameter \i 
can be adapted dynamically, yielding the Normalized 
LMS equation: 

60 Cjt + x =C k +aX k €{k)/{2{ZM+ 1)P} (Eq. 20) 

where 0<a<l. 

Pitch Tracking Methodology Using LMS Coefficient 
Computation 

The steps of the pitch tracking methodology were 
discussed above in a different order from the order in 
which steps are performed. Referring to FIG. 3, at 
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initialization, the process starts at step 200 by sampling 
the quasi-periodic input signal for N sampling periods 
(where NT, is larger than the period of the lowest pitch 
input signal in a predefined range of pitch values) and 
then using the AMDF technique to compute an initial 5 
period value To which is integer number P of sampling 
periods. Also, all the predictor coefficients are set to a 
value of zero, except c(0) which is set to a value of 1. 

After storing the new input signal sample for each 
successive sampling period (step 201), the following \q 
computations are performed, in sequence. First (at step 
202) the error signal €(k) is computed using the current 
signal sample values and the most recently computed 
predictor coefficients c(i). At step 204 the average sig- 
nal power x 2 is updated, and then at step 206 the predic- ^ 
tor coefficients c(i) are updated using equation 20 (see 
above). At step 208 the angular frequency o is com- 
puted based on the previously computed period value, 
and then the FIR filter's phase is computed at step 210 
using equation 5. At step 212, the FIR filter's phase 2 o 
delay is computed as the FIR filter's phase 0 divided by 
frequency <u. Finally, at step 214 the input signal's per- 
iod is computed by adding the computed phase delay to 
the integer time delay of the FIR filter P*T„ and the 
pitch is computed as the inverse of the computed per- 25 
iod. 

Steps 216 and 218 of the method are performed by the 
monitoring microprocessor in the preferred embodi- 
ment and concern the process of changing the integer 
delay period of the FIR filter as the pitch of the input 3Q 
signal increases or decreases. Whenever the integer 
delay period P of the FIR filter is changed by the digital 
signal processor (step 220), a "hold period" is started 
during which time further adjustment of the integer 
delay period P is prevented. The hold period would 35 
typically be between 50 and 100 sample periods for a 
system with 0.1 millisecond sample period (i.e., a 10 
KHz sampling clock). The hold period is invoked in 
order to prevent rapid oscillation of the integer delay 
period, by allowing the prediction coefficients to reach ^ 
their optimal values after each adjustment of the integer 
delay period P. Thus, at step 216, if a hold period is in 
effect, the process resumes at step 202 for processing the 
next sample period's data and time lagged data. 

If at step 216 it is determined that a hold period is not 
in effect, step 218 determines whether the magnitude of 
the FIR filter's phase delay is larger than a half sam- 
pling period (0.5T S ). If not, the process resumes running 
at step 202. If the phase delay's magnitude does exceed 
0.5T 5 , at step 220 the filter's integer delay period P is 
increased by one if the phase delay is larger than 0.5 
sample periods. Otherwise, at step 220 the filter's inte- 
ger delay period is less than —0.5 sample periods, and 
therefore the integer delay period is decreased by one. 
Then the prediction coefficients are swapped in accor- 
dance with equation 21, and a holding period is started, 
as discussed above. 

c0>±c(-i) for \i\ £M (Eq. 21) 

It should be noted that with a five element FIR filter, 60 
the predictor coefficients will accurately reflect phase 
delays of up about two sampling periods. To track even 
faster moving pitch changes, or to increase noise immu- 
nity, the number of FIR filter elements should be in- 
creased, for example, to seven. At step 220 the phase 65 
delay could be further tested to see if the magnitude of 
the phase delay is greater than 1.5 sampling periods, and 
if so, then the integer phase delay P would be increased 
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or decreased (in accordance with the sign of the phase 
delay) by a value of two. 

Another function performed by the monitoring mi- 
croprocessor is to monitor the pitch tracker's error 
signal for sudden increases. The error signal indicates 
when the pitch tracker has lost track of the input signal, 
such as when the pitch of the input signal suddenly 
changes by a large amount or otherwise changes faster 
than the pitch tracker can follow, or that the input 
signal has ceased to be quasi-periodic. When the error 
signal's value indicates that the pitch tracker has lost 
track of the input signal, the monitoring microprocessor 
sends commands to the DSP forcing it to restart is com- 
putation at step 200, where it recaptures the pitch of the 
input signal using the AMDF technique. 

To detect loss of tracking, the mean square error 
MSE (equal to the average of the square of the error 
signal €(k) over several sample periods) is compared 
with signal power x 2 . When error signal power (MSE) 
is greater than a preselected fraction of the signal power 
(i.e., MSE>/Jj?, where /3 is typically between 0.2 and 
0.8), the pitch tracker attempts to re-acquire a lock on 
the signal by returning to step 200. 

Alternate Pitch Tracking System Architectures 

Referring to FIG. 4, system 250 uses three parallel 
pitch trackers 252, 254 and 256 to avoid the settling time 
problems associated with the system shown in FIG. 2. 
Each of these pitch trackers is equivalent to the DSP 
implemented pitch tracker of FIG. 2, and the associated 
microprocessors are combined here into one monitoring 
microprocessor 260. The three pitch trackers are con- 
figured so that one pitch tracker 254 has an integer 
delay P0 while the other two pitch trackers have inte- 
ger delays of P0— 1 and P0+ 1 . All three pitch trackers 
252, 254 and 256 sample the input signal every period 
and compute the optimal prediction coefficients for that 
pitch tracker. As will be appreciated by those skilled in 
the art, some of the calculations for the three pitch 
trackers can be shared, allowing three pitch trackers to 
be implemented at relatively low computational cost. 

The monitoring microprocessor 260 monitors the 
FIR filter phase delay of center pitch tracker 254. If the 
phase delay is greater than 0.5T J} the integer delay in all 
three trackers is increased by one, the predictor coeffi- 
cients of the center pitch tracker with delay period P0 
(now increase by one) are copied to the P0— 1 pitch 
tracker 252. The predictor coefficients of the P0+1 
pitch tracker 256 are copied to the P0 pitch tracker 254, 
and the coefficients of the P0+ 1 pitch tracker 256 are 
time reversed as in equation 21 above. 

When the detected pitch of the input signal decreases, 
causing the phase delay of pitch tracker 254 to be less 
than — 0.5T„ a complimentary set of operations is per- 
formed: 

-Copy the coefficients of the P0 pitch tracker 254 to the 

P0+ 1 pitch tracker 256; 
-Copy the coefficients of the P0- 1 pitch tracker 252 to 

the PO pitch tracker 254; 
-Time reverse the coefficients of the P0— 1 pitch 

tracker 252. 

A hold period between coefficient swaps is required 
to avoid rapid oscillation of the delay line length when 
the input signal's period is near a 0.5 sample boundary, 
and also to allow the time reversed coefficients to adapt 
to optimal values. 



03/12/2004, EAST version: 1.4.1 



11 



5,353,372 



12 



In FIG. 5, there is shown a system 270 with an array 
of one hundred parallel pitch trackers 272 to 278 for 
monitoring a defined frequency range, such as 49 Hz to 
2500 Hz. The pitch trackers have preset integer delays 
varying in increments of two from a predefined low 5 
value (e.g., 4) associated with a maximum frequency 
(e.g., 2500 Hz) to a predefined high value (e.g., 202) 
associated with a minimum frequency (e,g., 49 Hz). A 
monitoring microprocessor 280 scans the error signals 
generated by all the pitch trackers, and selects the one 10 
with the smallest error signal and shortest period so as 
to select the pitch tracker most closely centered on the 
pitch of the input signal. The coefficients from the se- 
lected pitch tracker are read and converted into a phase 
delay, as described above, to derive an accurate repre- 
sentation of the pitch of the input signal. Once the fre- 
quency range of the system 270 is selected, the integer 
delays of the pitch trackers are not changed. Clearly, 
the number of pitch trackers needed in the array will 20 
depend on the frequency range being scanned. The 
number of pitch trackers used in the array could be 
reduced by changing the increment in delay periods 
between neighboring pitch trackers from 2T S to 3T* or 
«V 25 

In alternate embodiment, the invention could be im- 
plemented using an even number N of delayed signal 
samples and predictor coefficients. The use of an even 
number of predictor coefficients requires only minor 
adjustments to the process. For example, when tracking 30 
of the signal pitch requires changing the integer period 
P, equation 21 is applied by reversing the order of the 
predictor coefficients. 

While the present invention has been described with 
reference to a few specific embodiments, the descrip- 35 
tion is illustrative of the invention and is not to be con- 
strued as limiting the invention. Various modifications 
may occur to those skilled in the art without departing 
from the true spirit and scope of the invention as defined 
by the appended claims. 

What is claimed is: 

1. A pitch measurement system which determines the 
pitch of a quasi-periodic input signal, comprising: 

a signal sampler which samples said input signal's 
value x(t) during a sequence of sample periods; 

signal storage, coupled to said signal sampler, which 
stores a sequence of signal sample values represent- 
ing said input signal's value during said sequence of 
sample periods; ^ 

means, coupled to said signal storage, for generating 
an initial estimate P*T 5 of said input signal's period 
To where Tjis the length of each said sample period 
and P is an integer; 

means for storing a ordered set of N coefficient val- 55 
ues, where N is a positive integer; 

a finite impulse response (FIR) filter coupled to said 
signal storage for generating a predicted signal 
value x(k) for a sampling period k by multiplying 
each of N of said stored signal sample values, se- ^ 
lected in accordance with said initial estimate P*T* 
of said input signal's period, by corresponding ones 
of said N coefficient values and summing the result- 
ing values; 

eror means for generating an error signal e(k) corre- 65 
sponding to the difference between the actual 
stored signal value x(k) for sampling period k and 
the predicted signal value x(k); 
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coefficient updating means for periodically updating 
said coefficient values so as to minimize said error 
signal e(k)'s mean square value; and 

measurement output means for generating an output 
measurement equal to P-T^plus a phase delay value 
that is equal to a predefined function of said up- 
dated coefficient values. 

2. The pitch measurement system of claim 1, wherein 
said coefficient updating means periodically updates 
each of said N coefficient values by an amount propor- 
tional to the error signal e(k) multiplied by a corre- 
sponding one of the N selected signal values: 

Cjk + i=Cjfc+2pjr*c(A) 

where C* is a vector containing the N coefficient values 
for sampling period k, 2\i is an adaption parameter that 
controls the rate at which the coefficient values are 
adjusted, and X* is a vector containing the N selected 
stored signal sample values. 

3. The pitch measurement system of claim 1, wherein 
said measurement output means generates said phase 
delay value for each sampling period k with respect to 
a previously computed output measurement. 

4. The pitch measurement system of claim 1, further 
including tracking means for increasing the value of P 
when said phase delay value is larger than 0.5T A for 
decreasing the value of P when said phase delay value is 
less than — 0.5Tj> and for swapping said coefficient 
values so as to reverse said ordered set of N coefficient 
values whenever the value of P is changed. 

5. A pitch measurement system which determines the 
pitch of a quasi-periodic input signal, comprising: 

a signal sampler which samples said input signal's 
value x(t) during a sequence of sample periods; 

signal storage, coupled to said signal sampler, which 
stores a sequence of signal sample values represent- 
ing said input signal's value during said sequence of 
sample periods; 

means, coupled to said signal storage, for generating 
an initial estimate P»T,of said input signal's period 
To where T^is the length of each said sample period 
and P is an integer; 

means for storing a set of 2M+1 coefficient values 
c(i), for i— — M to -fM, where M is a positive 
integer; 

a finite impulse response (FIR) filter coupled to said 
signal storage for generating a predicted signal 
value x(k) for sampling period k in accordance 
with the following formula: 



where x(k— P+0 is the stored signal sample value for 

sampling period k— F+i; 
error means for generating an error signal «(k) corre- 
sponding to the difference between the actual 
stored signal value x(k) and the predicted signal 
value x(k); 

coefficient updating means for periodically updating 
said coefficient values so as to minimize said error 
signal €(k)'s mean square value; and 

measurement output means for generating an output 
measurement equal to P-T 5 plus a phase delay value 
that is equal to a predefined function of said up- 
dated coefficient values. 
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6. The pitch measurement system of claim 5, wherein 
said coefficient updating means periodically updates 
said coefficient values by an amount proportional to the 
error signal e(k) multiplied by each of the 2M+ 1 se- 
lected stored signal sample values: 5 

c* + i=c*+2|aw*) 

where C* is a vector containing the 2M+ 1 coefficient 
values for sampling period k, 2^ is an adaption parame- io 
ter that controls the rate at which the coefficient values 
are adjusted, and X* is a vector containing the 2M+ 1 
selected stored signal sample values for sample periods 
k-P-Mtok-P+M 

7. The pitch measurement system of claim 5, wherein IS 
said measurement output means generates said phase 
delay value for each sampling period k with respect to 

a previously generated output measurement 

8. The pitch measurement system of claim 5, further 
including tracking means for increasing the value of P 20 
when said phase delay value is larger than Q.5T S , for 
decreasing the value of P when said phase delay value is 
less than — O.STj, and for swapping each of said coeffici- 
ent values c(i) with coefficient value c(— i) whenever 
the value of P is changed. 25 

9. A pitch measurement system which determines the 
pitch of a quasi-periodic input signal, comprising: 

a signal sampler which samples said input signal's 
value x(t) during a sequence of sample periods; 

signal storage, coupled to said signal sampler, which 30 
stores a sequence of signal sample values represent- 
ing said input signal's value during said sequence of 
sample periods; 

a multiplicity of parallel signal pitch trackers, each 
said signal pitch tracker having assigned thereto an 35 
initial estimate P-T* of said input signal's period To 
where T s is the length of each said sample period 
and P is an integer; wherein a different value of P 
is assigned to each of said multiplicity of parallel 
signal pitch trackers; each said signal pitch tracker 40 
generating an estimated signal period value and an 
error signal value; 

monitoring means, coupled to said multiplicity of 
parallel signal pitch trackers, for monitoring said 
error signal values generated by said multiplicity of 45 
parallel signal pitch trackers and for outputting the 
estimated signal period value generated by the one 
of said signal pitch trackers with the lowest error 
signal value; 

each of said multiplicity of parallel signal pitch track- 50 
ers including: 

means for storing a ordered set of N coefficient 
values, where N is a positive integer; 

a finite impulse response (FIR) filter coupled to 
said signal storage for generating a predicted 55 
signal value x(k) for a sampling period k by mul- 
tiplying each of N of said stored signal sample 
values, selected in accordance with said initial 
estimate P-T, of said input signal's period, by 
corresponding ones of said N coefficient values 60 
and summing the resulting values; 

error means for generating an error signal c(k) 
corresponding to the difference between the 
actual stored signal value x(k) for sampling per- 
iod k and the predicted signal value x(k); 65 

coefficient updating means for periodically updat- 
ing said coefficient values so as to minimize said 
error signal e(k)'s mean square value; and 
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output means for generating an output measurement 
equal to P Tj plus a phase delay value that is equal 
to a predefined function of said updated coefficient 
values. 

10. A pitch measurement method which determines 
the pitch of a quasi-periodic input signal, the steps of the 
method comprising: 

sampling said input signal's value x(t) during a se- 
quence of sample periods; 

storing a sequence of signal sample values represent- 
ing said input signal's value during said sequence of 
sample periods; 

using said stored sequence of signal sample values, 
generating an initial estimate P-T* of said input 
signal's period To where Tj is the length of each 
said sample period and P is an integer; 

storing an ordered set of N coefficient values, where 
N is a positive integer; 

selecting, in accordance with said initial estimate 
P-Tj of said input signal's period, N of said stored 
signal sample values; (B) generating a predicted 
signal value x(k) for a sampling period k by multi- 
plying each of said N selected signal sample values 
by corresponding ones of said N coefficient values 
and summing the resulting values; and (C) generat- 
ing an error signal €(k) corresponding to the differ- 
ence between the actual stored signal value x(k) for 
sampling period k and the predicted signal value 
x(k); 

updating said coefficient values so as to minimize said 
error signal e(k)'s means square value; and 

generating an output measurement equal to P-Tjplus 
a phase delay value that is equal to a predefined 
function of said updated coefficient values. 

11. The pitch measurement method of claim 10, 
wherein said coefficient value updating step periodi- 
cally updates each of said coefficient values by an 
amount proportional to the error € multiplied by a cor- 
responding one of the N selected signal values: 

C k+ i = C k +2itXk€(k) 

where C* is a vector containing the N coefficient values 
for sampling period k, 2fi is an adaption parameter that 
controls the rate at which the coefficient values are 
adjusted, and X* is a vector containing the N selected 
stored signal sample values. 

12. The pitch measurement method of claim 10, 
wherein said output measurement generating step com- 
putes for each sampling period k said phase delay value 
with respect to a previously computed output measure- 
ment. 

13. The pitch measurement method of claim 10, fur- 
ther including the steps of: 

increasing the value of P when said phase delay value 
is larger than 0.5T S ; decreasing the value of P when 
said phase delay value is less than — 0.5T 5 ; and for 
swapping said coefficient values so as to reverse 
said ordered set of N coefficient values whenever 
the value of P is changed. 

14. A pitch measurement method which determines 
the pitch of a quasi-periodic input signal, the steps of the 
method comprising: 

sampling said input signal's value x(t) during a se- 
quence of sample periods; 

storing a sequence of signal sample values represent- 
ing said input signal's value during said sequence of 
sample periods; 



03/12/2004, EAST Version: 1.4.1 



15 



5,353,372 



16 



using said stored sequence of signal sample values, 
generating an initial estimate P-T* of said input 
signal's period To where T 5 is the length of each 
said sample period and P is an integer; 

storing a set of 2M+ 1 coefficient values c(i), for 
i = — M to -f M, where M is a positive integer; 

selecting 2M+1 of said stored signal sample values 
corresponding to sample periods delayed from a 
current signal sample period k by periods ranging 
from (P— M)T, to (P+M)T* generating a pre- 
dicted signal value x(k) for sampling period k in 
accordance with the following formula: 



. 2 rf* - * + i*<0 
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where x(k— P+i) is the stored signal sample value for 
sampling period k— P+i, and generating an error signal 
c(k) corresponding to the difference between the actual 20 
stored signal value x(k) and the predicted signal value 
x(k); 

updating said coefficient values so as to minimize said 
error signal e(k)'s means square value; and 

generating an output measurement equal to P-T 5 plus 
a phase delay value that is equal to a predefined 
function of said updated coefficient values. 
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15. The pitch measurement method of claim 14, 
wherein said coefficient value updating step periodi- 
cally updates said coefficient values by an amount pro- 
portional to the error e multiplied by each of the 2M+ 1 
time lagged signal values: 

where c* is a vector containing the 2M+ 1 coefficient 
values for time period k, 2p, is an adaption parameter 
that controls the rate at which the prediction coeffici- 
ents are adjusted, and Xjt is a vector containing the 
2M-f 1 selected stored signal sample values for sample 
periods k— P-M to k— P+M 

16. The pitch measurement method of claim 14, 
wherein said output measurement generating step com- 
putes for each sampling period k said phase delay value 
with respect to a previously computed estimate signal 
period value for said input signal. 

17. The pitch measurement method of claim 14, fur- 
ther including the steps of: 

increasing the value of P when said phase delay value 
is larger than 0.5T 5 ; decreasing the value of P when 
said phase delay value is less than — 0.5T r ; and 
swapping each of said coefficient values c(i) with 
coefficient value c(— i) whenever the value of P is 
changed. 
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